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Wo measure the mass dependence of the forward-backward charge asymmetry in 157,553 pp — > 
Z/7* — > e+e~ interactions, corresponding to 5.0 fb"'^ of integrated luminosity collected by the 
DO experiment at the Fermilab Tevatron Collider at \fs = 1.96 TeV. The effective weak mixing 
angle (Oig) from this process involving predominantly the first generation of quarks is extracted 
as sin'^eis = 0.2309 ± 0.0008 (stat.) ± 0.0006 (syst.). We also present the most precise direct 
measurement of the vector and axial-vector couplings of u and d quarks to the Z boson. 

PACS numbers: 12.15.Ji, 12.15.Mn, 13.85.Qk, 14.70.Hp, 13.38.Dg 



INTRODUCTION 

Electron-positron pairs (e+e~) can be produced 
through the Drell-Yan process over a large invariant mass 
range at the Fermilab Tevatron Collider. In the standard 
model (SM) of particle physics, the process occurs to first 
order via qq annihilation into a real (or virtual) Z boson 
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or a virtual photon (7*). While the coupling of a fermion 
(/) to the photon is purely a vector coupling, the cou- 
pling of the same fermion to the Z boson has both vector 
{gy) and axial-vector (g^) components: 

9v = li -2qfsm^ Ow, /-.x 

„/ _ rf 
9 a — ^3^ 

£ 

whore and qj arc the third component of the weak 
isospin and the charge of the fermion, and 9w is the weak 
mixing angle [1]. The presence of both vector and axial- 
vector couplings gives rise to an asymmetry in the dis- 
tribution of the polar angle 0* of the negatively charged 
lepton relative to the incoming quark direction in the rest 
frame of the lepton pair. To minimize the effect of the 
unknown transverse momentum of the incoming quarks, 
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we calculate 6* in the Collins-Soper reference frame [2] 
as 



(Pi-'Pf-Pi-PT^)^ (2) 



\QWQ' + Qt ' ' ' ' 

where Q {Qt) is the four momentum (transverse momen- 
tum) of the lepton pair, and Pi and Pj are the four mo- 
menta of the lepton and anti-lepton, respectively. They 
are measured in the lab frame, and the momenta P^ are 
defined as 



1 



(3) 



where P° and Pj^ are the energy and the longitudinal 
component of the lepton momentum, respectively. In 
the Collins-Soper frame, the polar axis is defined as the 
bisector of the proton beam momentum Pi and the nega- 
tive of the anti-proton beam momentum, —P2, when the 
proton and anti-proton are boosted into the rest frame 
of the lepton pair, as shown in Fig. 1 [3]. 
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FIG. 1: The Collins-Soper reference frame. The bisector of 
the proton beam momentum P\ and the negative of the anti- 
proton beam momentum — P2 are used to measure the angle 
d* . The momenta P\ and P2 are measured in the e^e~ rest 
frame. 

Events with electron cos 0* > are classified as forward 
(F), and those with electron cos0* < are classified as 
backward (P). The forward-backward charge asymme- 
try, AfBi is defined by 



ApB = 



ap + ctb' 



(4) 



where ap and aB are the cross sections for forward and 
backward processes, respectively. 

The SM leading order (LO) prediction [4] for ApB as a 
function of the dielectron invariant mass {M^e) is shown 
in Fig. 2 for uu -> Z/j* — > e+e~, d J — > Z/j* — > e+e~, 
and pp — > Z/j* — > e+e^ with the CTEQ6L1 parton 
distribution functions (PDFs) [5]. Around the Z pole, 
the asymmetry is proportional to both the vector and 
axial-vector couplings of the Z boson to the fermions 
and is numerically close to 0. At large invariant mass, 
the asymmetry is dominated by Z/j* interference and is 
almost constant (sa 0.6). In the high mass region, the 



ApB measurement can be used to investigate possible 
new phenomena that may alter ApB, such as new neutral 
gauge bosons or large extra dimensions [6-14]. 
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FIG. 2: [color online] The SM LO Afb prediction as a 
function of the dielectron invariant mass for uu — >■ e'^e', 
dd e'^ e~ , and pp ^ e~ [4]. 

In the vicinity of the Z pole, ApB is sensitive to the 
effective weak mixing angle (sin^ 9^^g) for each ferniion 
species, /, involved in a particular measurement. To all 
orders in perturbation theory [1, 15], sin^ is related to 
the vector and axial-vector couplings by the expression 



9v/9a 



1 - 4|g/| sin^ 



(5) 



This charged lepton effective mixing angle sin^ 

varies as a function of the momentum transfer at which 

it is measured. Conventionally, it is quoted at the Z pole 

[sin^ 6^^ff{Mz)] , and it is identical for e, ^, and r leptons, 

due to lepton universality. 

In the SM, asymmetries measured at the Z pole [15] 

2 f 

depend only on the value of sin 9^ff for the fermions be- 
ing considered. Because of the small ratio of vector and 
axial-vector couplings for leptons, the sensitivity of lep- 
tonic asymmetries to the changes in effective mixing an- 
gle arises predominantly through the variation of the lep- 
tonic couplings and not those of the quarks. Therefore, 
it is customary to express ApB measurements in terms 
of sin 9^g. In order to extract sin from ApB under 
a consistent SM definition and compare results with pre- 
vious measurements, we take into account the difference 
between the electroweak radiative corrections for elec- 
trons and u/d quarks using the relations [15-17] 



sin^ 6i^ff = sin2 61^^- 0.0001, 



'^cft = sin2 0. 



off 



0.0002. 



(6) 



Precise determinations of sin have been made in 
many processes at different scales. They include 
atomic parity violation {\Q^\ ~ 10" GeV^) [18], M0ller 
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scattering using a polarized electron beam and unpolar- 
ized target {\Q^\ 0.03 GeV^) [19], the NuTeV deep 
inelastic neutrino and anti-neutrino scattering on iron 
(IQ^I « 4 GcV^) [20], and a number of measurements 
employing e+e" collisions by the LEP and SLD Col- 
laborations (IQ^I Ri M|) [15]. The current world av- 
erage value of sin^ ff^g is 0.23153±0.00016 [15]. The two 
most precise determinations of sin^ come from the b- 
quark forward-backward asymmetry at LEP, vl^^ , with 
sin^ 0.23221 ± 0.00029, and the left-right asymme- 

try at SLD, AiriSLB), with sin^ e^g = 0.23098 ±0.00026. 
These two measurements differ from each other by about 
three standard deviations, and deviate by +2.1 standard 
deviations and —1.8 standard deviations from the global 
fit, respectively. 

The LEP Collaborations also measured sin^ 6*^^. from 
the inclusive hadronic charge asymmetry {Q^}?g)- with 
larger uncertainties governed by the ambiguity of charge 
separation for final state quark species. Furthermore, 
the hadronic charge asymmetry arising from u- and d- 
type quarks are in opposite directions, partially cancel- 
ing. Thus, modifications to the SM that would affect only 
u and d quark couplings are poorly constrained. Drell- 
Yan processes at hadron colliders, in which the initial 
state is dominated by the light u and d quarks in the pro- 
ton, provide a much less ambiguous measurement of the 
light quark couplings. The dominant systematic uncer- 
tainty at the Tevatron originates from the quark compo- 
sition of the proton and anti-proton, which has been well 
constrained and parametrized by PDFs [5]. The use of 
the CoUins-Soper frame reduces possible confounding ef- 
fects from higher order quantum chromo dynamics (QCD) 
corrections. 

Previous direct measurements of u and d quark cou- 
plings to the Z boson arc of limited precision ([15, 21, 
22]). With precise determination of the leptonic cou- 
plings from LEP and SLD, we can interpret the mea- 
surement of the forward-backward asymmetry directly 
in terms of the vector and axial- vector couplings of the 
u/d quarks. 

At the Tevatron, measurements of the Afb, sin^ 0^^, 
gy and g^"^^ have been performed by the CDF and DO 
Collaborations [21, 23-25]. The largest integrated lumi- 
nosity used for these measurements was 1.1 fb~^ for Afb 
and sin^ 6^^ measurements [25], and 72 pb~^ for ''^ and 
g^"^^ measurements [21]. In this Article we present new 

measurements of the quantities sin^ gy "^^ and g^'^'' 
based on 5.0 fb~^ of integrated luminosity [26], collected 
using the DO detector [27] between April 2002 and April 
2009. 



APPARATUS AND EVENT SELECTION 

The DO detector [27] comprises a central tracking sys- 
tem, a calorimeter and a muon system. The central track- 
ing system is composed of a silicon microstrip tracker 
(SMT) and a central fiber tracker (CFT), both located 
within a 2 T superconducting solenoidal magnet and op- 
timized for tracking and vertexing capabilities at detector 
pseudorapidities of |77det| < 3 [28]. Three liquid argon and 
uranium calorimeters provide coverage of |r?det| < 3.2 for 
electrons with gaps between cryostats creating an ineffi- 
cient electron detection region between 1.0 < |77det| < 1.5. 
The electromagnetic (EM) section of the calorimeter is 
segmented into four longitudinal layers with transverse 
segmentation of A?; x A(j) = 0.1 x 0.1, except for the 
third layer, where it is 0.05 x 0.05. The calorimeter is well 
suited for a precise measurement of electron and photon 
energies, providing a resolution of ~ 3.6% at an incident 
energy of ~ 50 GeV. The muon system surrounds the 
calorimctry and consists of three layers of scintillators 
and drift tubes and 1.8 T iron toroids with a coverage of 
jr^detl < 2. The three-level trigger system and data ac- 
quisition systems arc designed to accommodate the high 
instantaneous luminosity of Run II. A logical OR of di- 
electron triggers is used to collect the data, resulting in 
a trigger efficiency close to 100% for signal events that 
passed the offline selection. 

To select Zf-y* — > e+e" events, we require two EM 
shower candidates with transverse energy Et > 25 GeV 
measured in the calorimeter. An isolation cut is im- 
posed on the candidates, requiring that the fraction of 
their energy in an annular central (endcap) calorime- 
ter cone of radius 0.2 < ATZ < 0.4 must be less than 
15% (10%) of the energy in the cone of A7^ < 0.2, 
where ATZ = (Atj)^ -|- (A^)^. The candidates are fur- 
ther required to have a significant fraction of their en- 
ergy deposited in the EM calorimeter compared to that 
in the hadron calorimeter, and to have a shower shape 
consistent with that expected for an electron. At least 
one electron candidate is required to be in the central 
(l%et| < 1-0) fiducial region and spatially matched to 
a reconstriicted track, while the other candidate may 
be either in the central or endcap (1.5 < |r?det| < 2.5) 
calorimeter. No track requirement is imposed on can- 
didates in the endcap calorimeter, since the track re- 
construction efficiency is degraded in this region. If an 
event has both candidates in the central calorimeter (CC 
events), the two candidates are further required to have 
opposite charges. For events with one candidate in the 
central and the other in the endcap calorimeter (CE 
events), the charge of the central EM candidate is used 
to determine if it is a forward or a backward event. To 
suppress multijet background in CE events, the electron 
candidates in the endcap calorimeter are required to pass 
isolation criteria in the tracker, requiring the scalar sum 
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of the transverse momenta of tracks in the annulus 0.05 
< ATZ < 0.4 centered around the electron direction to be 
smaUer than 1.5 GeV. Events are further required to have 
the reconstructed pp interaction vertex within 40 cm of 
the detector center in the coordinate along the z axis and 
a reconstructed invariant mass of the electron pair {Mee) 
between 50 and 1000 GeV. 

A total of 157,553 events remain after application of 
all selection criteria, with 73,755 CC events and 83,798 
CE events. The forward-backward charge asymmetries 
are measured in 15 M^e bins in the range 50 < M^e < 
1000 GeV. The bin widths are chosen considering the 
statistics of the sample and the mass resolution of the 
detector. The bin widths and the numbers of forward and 
backward events for each mass bin are listed in Table I. 



SIGNAL AND BACKGROUNDS 
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FIG. 3: [color online] The data/MC ratio of electron selection 
efficiencies as a function of invariant mass, for forward and 
backward events. 



Monte Carlo (MC) samples for the Z/7* -J> e+e" pro- 
cess are generated using the PYTHIA event generator with 
CTEQ6L1 PDFs, followed by a detailed GEANT-based 
simulation [29] of the DO detector response. This simu- 
lation is then improved by corrections for observed defi- 
ciencies in the detector simulation and for higher order 
physical effects not included in pythia. 

The event selection criteria can result in different effi- 
ciencies for forward and backward events. The electron 
selection efficiencies are independently measured from 
Zf-f* e+e~ events in data and in the MC, where one 
electron is selected in the central calorimeter using tight 
calorimeter shower shape cuts and track quality cuts, and 
the second electron is used as a probe to determine the 
detection efficiencies. These efficiencies are measured for 
forward and backward events separately. For data, the 
background in each mass bin is estimated and subtracted 
prior to the measurement of the efficiencies. The ratios 
between data and MC electron selection efficiencies for 
forward and backward candidates as a function of Mee for 
electrons in the central calorimeter are shown in Fig. 3. 
The ratios are constant within statistical uncertainties, 
with the largest deviations observed in a few mass bins 
around 70 and 130 GeV. Efficiency corrections derived 
using data presented on Fig. 3 are applied to the MC sep- 
arately for forward and backward events to account for 
the mis-modeling of electrons' shower shapes and track 
matching efficiencies. In addition, the MC is adjusted to 
reproduce the calorimeter energy scale and resolution, as 
well as the distributions of the instantaneous luminosity 
and the event vertex position observed in data. 

Next-to-next-to-leading order (NNLO) QCD correc- 
tions [30] for Z/j* boson production are applied to the 
simulated pythia sample by reweighting the Mee dis- 
tribution, and non-perturbative and next-to-leading or- 
der (NLO) corrections by reweighting the Z/j* boson 
transverse momentum and rapidity distributions [17, 31]. 



The effective mixing angle must be corrected to include 
higher order quantum electrodynamics (QED) and weak 
interaction contributions that are not present in our MC. 
These higher order corrections are determined using the 
ZGRAD2 program [16]. 

The largest background originates from multijet events 
in which jets are mis-reconstructed as electrons. Smaller 
background contributions arise from other SM processes 
that produce at least one real electron or photon in 
the final state. SM backgrounds, such as Z/j* — >■ tt, 
W + X, WW, WZ, 77, and ti, are estimated using the 
MC. Higher order corrections to the cross sections have 
been applied [31-33]. The multijet background is esti- 
mated using collider data by fitting the Mee distribution 
in the Z pole region (with other SM backgrounds sub- 
tracted) to the sum of the shape predicted by the cor- 
rected Z/j* — >■ e+e" signal MC and the shape measured 
from a multijet-enriched sample. The multijet-enriched 
sample is selected by reversing the shower shape require- 
ment on the two electron candidates. The average mul- 
tijet background fraction over the entire mass region is 
found to be approximately 0.9%. The numbers of back- 
ground events from each source in forward and backward 
samples are listed in Tables II and III. 

Comparisons of data and the sum of signal and back- 
ground for Mee and cos 6** are shown in Figs. 4 and 5. In 
the Mee bin 450-500 GeV, the data differ from the SM 
prediction by 1.8 standard deviations. Reasonable agree- 
ment is observed for all distributions in both forward and 
backward samples for all 15 Mee bins. The CC and CE 
raw Afb (not yet unfolded) distributions as functions 
of Mee are then calculated from background-subtracted 
data. 



7 



TABLE I: Numbers of forward and backward CC and CE events in each Mee bin after all selections. 



Forward Backward Forward Backward 



50 — 


60 


276 


319 


54 


70 


60 — 


70 


464 


711 


238 


413 


70 — 


75 


411 


545 


285 


495 


75 — 


81 


852 


1062 


778 


1240 


81 — 


86.5 


3359 


3559 


3804 


4245 


86.5- 


89.5 


6681 


6642 


8339 


7591 


89.5- 


92 


9297 


8717 


11098 


9534 


92- 


97 


12076 


11109 


14281 


11412 


97- 


105 


2890 


2173 


3711 


2150 


105- 


115 


680 


431 


1125 


395 


115- 


130 


408 


189 


764 


229 


130- 


180 


439 


150 


845 


269 


180- 


250 


138 


61 


229 


73 


250- 


500 


63 


45 


86 


24 


500- 


1000 


7 


1 


1 






TABLE IL Estimated number of background events in each Mee bin in the forward sample. 
Mee (GeV) Z/-y* tt W + X WW WZ 77 tt Multijit 



50- 


60 


12.5 


± 


0.91 


11.7 


± 


4.41 


1.65 


± 


0.14 


0.14 


± 


0.01 


1.02 


± 


0.35 


1.00 


± 


0.11 


38.6 


± 


0.04 


60- 


70 


29.4 


± 


1.44 


20.4 


± 


7.17 


3.29 


± 


0.24 


0.22 


± 


0.01 


3.10 


± 


0.43 


1.83 


± 


0.21 


105. 


± 


0.10 


70- 


75 


16.6 


± 


0.97 


17.2 


± 


4.38 


1.68 


± 


0.14 


0.18 


± 


0.01 


1.08 


± 


0.26 


1.09 


± 


0.13 


69.2 


± 


0.09 


75- 


81 


14.5 


± 


0.91 


16.6 


± 


4.86 


1.55 


± 


0.13 


0.31 


± 


0.01 


1.59 


± 


0.28 


1.37 


± 


0.15 


85.8 


± 


0.10 


81- 


86.5 


5.16 


± 


0.72 


21.4 


± 


8.33 


1.82 


± 


0.14 


0.80 


± 


0.03 


2.26 


± 


0.31 


1.16 


± 


0.14 


80.3 


± 


0.10 


86.5- 


89.5 


0.94 


± 


0.49 


8.03 


± 


2.67 


1.10 


± 


0.12 


1.72 


± 


0.06 


1.16 


± 


0.26 


0.56 


± 


0.07 


40.6 


± 


0.07 


89.5- 


92 


1.63 


± 


0.69 


9.73 


± 


3.07 


0.86 


± 


0.11 


2.79 


± 


0.11 


0.29 


± 


0.24 


0.68 


± 


0.08 


31.0 


± 


0.07 


92- 


97 


1.04 


± 


0.49 


18.9 


± 


5.34 


1.64 


± 


0.14 


3.66 


± 


0.13 


2.00 


± 


0.30 


0.99 


± 


0.12 


62.7 


± 


0.10 


97- 


105 


1.38 


± 


0.14 


24.4 


± 


10.8 


2.72 


± 


0.17 


0.80 


± 


0.03 


2.00 


± 


0.31 


1.50 


± 


0.15 


88.5 


± 


0.12 


105- 


115 


1.21 


± 


0.11 


23.4 


± 


11.3 


3.03 


± 


0.19 


0.33 


± 


0.01 


1.91 


± 


0.30 


1.22 


± 


0.12 


94.5 


± 


0.12 


115- 


130 


1.33 


± 


0.11 


30.0 


± 


15.5 


4.33 


± 


0.25 


0.38 


± 


0.01 


2.00 


± 


0.29 


2.16 


± 


0.19 


108. 


± 


0.13 


130- 


180 


2.38 


± 


0.50 


41.3 


± 


27.6 


9.87 


± 


0.55 


0.82 


± 


0.03 


4.51 


± 


0.27 


3.94 


± 


0.40 


174. 


± 


0.17 


180- 


250 


0.57 


± 


0.03 


17.1 


± 


14.5 


4.53 


± 


0.26 


0.48 


± 


0.02 


2.84 


± 


0.18 


1.73 


± 


0.17 


42.5 


± 


0.08 


250- 


500 


0.19 


± 


0.02 


4.80 


± 


3.39 


1.93 


± 


0.15 


0.25 


± 


0.01 


1.24 


± 


0.12 


0.53 


± 


0.06 


8.70 


± 


0.04 


500- 


1000 




< 


0.01 




< 


0.01 


0.07 


± 


0.00 




< 


0.01 


0.04 


± 


0.03 




< 


0.01 


0.02 


± 


0.00 



TABLE 111: Estimated rmrribcr of background events in each Mee bin in the backward sample. 



Mee 


(GeV) 


Z/j 




-> TT 


w + x 


WW 


WZ 




77 






tt 




Multijct 


50 


-60 


7.52 


± 


0.80 


4.27 


± 


1.83 


2.90 


± 


0.20 


0.20 


± 


0.01 


0.54 


± 


0.25 


1.40 


± 


0.16 


38.8 


± 


0.04 


60 


-70 


26.2 


± 


1.36 


4.20 


± 


2.39 


3.73 


± 


0.24 


0.28 


± 


0.01 


2.85 


± 


0.35 


2.06 


± 


0.24 


108. 


± 


0.10 


70 


-75 


12.4 


± 


0.85 


3.06 


± 


1.33 


2.64 


± 


0.18 


0.18 


± 


0.01 


2.32 


± 


0.32 


1.12 


± 


0.12 


70.2 


± 


0.09 


75 


-81 


5.13 


± 


0.72 


2.74 


± 


1.30 


2.66 


± 


0.18 


0.34 


± 


0.01 


1.66 


± 


0.28 


1.21 


± 


0.14 


85.9 


± 


0.10 


81 


-86.5 


1.10 


± 


0.50 


5.30 


± 


2.49 


1.68 


± 


0.13 


0.76 


± 


0.03 


2.33 


± 


0.32 


1.08 


± 


0.13 


78.8 


± 


0.11 


86.5 


-89.5 


0.53 


± 


0.49 


4.25 


± 


1.56 


1.66 


± 


0.14 


1.59 


± 


0.06 


1.44 


± 


0.27 


0.45 


± 


0.07 


39.5 


± 


0.08 


89.5 


-92 


0.15 


± 


0.09 


3.80 


± 


1.67 


1.32 


± 


0.13 


2.43 


± 


0.09 


1.09 


± 


0.26 


0.53 


± 


0.07 


31.5 


± 


0.07 


92 


-97 


0.24 


± 


0.10 


2.34 


± 


0.94 


2.04 


± 


0.14 


3.17 


± 


0.11 


2.68 


± 


0.34 


1.05 


± 


0.12 


64.1 


± 


0.10 


97 


-105 


0.30 


± 


0.03 


6.98 


± 


2.64 


3.93 


± 


0.24 


0.75 


± 


0.02 


2.20 


± 


0.31 


1.74 


± 


0.20 


90.0 


± 


0.12 


105 


-115 


0.26 


± 


0.03 


5.47 


± 


2.63 


3.07 


± 


0.19 


0.35 


± 


0.01 


2.72 


± 


0.34 


1.53 


± 


0.15 


97.3 


± 


0.12 


115 


-130 


0.46 


± 


0.10 


6.22 


± 


3.89 


3.25 


± 


0.20 


0.37 


± 


0.01 


2.60 


± 


0.33 


2.30 


± 


0.24 


110. 


± 


0.13 


130 


-180 


0.76 


± 


0.49 


17.0 


± 


10.2 


5.58 


± 


0.29 


0.63 


± 


0.03 


4.97 


± 


0.39 


4.08 


± 


0.40 


170. 


± 


0.17 


180 


-250 


0.30 


± 


0.48 


3.49 


± 


2.92 


2.42 


± 


0.16 


0.28 


± 


0.01 


3.01 


± 


0.18 


1.70 


± 


0.17 


46.3 


± 


0.08 


250 


- 500 


0.04 


± 


0.01 


0.62 


± 


0.85 


0.58 


± 


0.11 


0.08 


± 


0.00 


1.28 


± 


0.20 


0.50 


± 


0.06 


8.90 


± 


0.04 


500 


- 1000 




< 


O.Oi 




< 


0.01 




< 


0.01 




< 


0.01 


0.0:i 




0.00 




< 


0.01 


0.02 




0.00 
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FIG. 4: [color online] Comparisons of the dielectron invari- 
ant mass between data and the sum of signal and background 
predictions for combined CC and CE events. The insert fo- 
cuses on the Z pole region from 50 GeV to 130 GeV, where 
good agreement between data and the sum of signal and back- 
ground predictions is essential to perform the unfolding. 




FIG. 5: [color online] Comparisons of the cosS* between data 
and the sum of signal and background predictions for com- 
bined CC and CE events. 

DETECTOR RESOLUTION AND ACCEPTANCE 

The finite energy resolution in the determination of the 
track curvature may result in the assignment of events in 
different bins of invariant mass and in changes in the 
forward/backward classifications. These bin migration 
effects in the raw ApB distribution are corrected using 
an unfolding procedure based on the iterative application 
of the matrix inversion method [34] , as in a previous DO 
analysis [25] . The CC and CE raw ApB distributions are 
unfolded separately and then combined. We correct for 
both the wrong classification in terms of dielectron invari- 
ant mass and for the wrong forward/backward assign- 
ment by defining four detector response matrices. The 



response matrix represents the probability that an 
event which at the generator level was classified to corre- 
spond to forward scattering angles and to have Mge in the 
j-th bin, to be reconstructed in the z-th bin in Mee, with- 
out any change to the forward/backward assignment. We 
similarly define matrix for events which are classi- 
fied as backward both at the generator and at the recon- 
struction level, and the i?,^^ and Rff matrices for events 
in which the forward/backward assignment changes due 
to detector and reconstruction effects. Tests of the un- 
folding procedure are performed comparing the generator 
level distribution with the one obtained after unfolding 
the events processed through the full detector simulation 
and reconstructed as data. 

The bin purity is defined as the ratio between the num- 
ber of events generated in a mass bin and also recon- 
structed in the same mass bin (A^gg™) ^'^d the number of 
events reconstructed in this mass bin (TV"™). The low- 
est purity occurs for the two mass bins below the Z pole 
(81 < Mee < 86.5 GeV and 86.5 < Mee < 89.5 GeV) and 
is about 25%. Since the corrected MC can describe the 
data mass spectra with finer binning as shown in Fig. 4, 
these low purity bins can be well modeled. The rest of 
the mass bins have purity varying between 50% and 96%. 

After unfolding for detector resolution effects, the data 
are further corrected for acceptance. Using the corrected 
signal MC, we derive corrections for kinematic and geo- 
metric acceptance and for electron identification efficien- 
cies. 



CHARGE MISIDENTIFICATION RATE 

The electron charge determines if an event is forward 
or backward. Mismeasurement of the sign may result in a 
dilution oi Aps- The charge misidentification probability 
Jq is given by 

/, = ^^s5/(iVss + A^os), (7) 

where N^s (Nqs) is the total number of Z/'j* — e+e~ 
events reconstructed with same-sign (opposite-sign) elec- 
trons. Since few same-sign events are observed in data, 
the corrected Z/j* MC is used to determine the shape of 
the misidentification probability as a function of dielec- 
tron invariant mass. The overall normalization is set by 
the misidentification probability determined from data 
around the Z pole. The misidentification probability is a 
function of Mee and rises from 0.75% at Mee = 50 GeV to 
3.2% for Mee > 500 GeV. The charge misidentification 
probability is included as a dilution factor T) in Afb, 
with V = {1 ~ 2fq)/{l - 2fq + /2) for CC events and 
V = {1- 2fg) for CE events. 
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SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties in the measurement of 
sin^^gjj, and of the unfolded distribution of ApB, gv, 
and qa, are Usted below. 

• PDF 

Uncertainties in the input parton distributions lead 
to uncertainties in the event acceptance. The sys- 
tematic uncertainty due to the PDFs uncertainty 
is estimated by reweighting the central PDFs us- 
ing the 40 CTEQ6.1M error sets, and the 90% C.L. 
uncertainty is calculated using the prescription sug- 
gested by the CTEQ group [5]. 

• Electron energy scale and resolution 

The energy scale and resolution for electrons in 
MC are corrected to match the observed Z boson 
pole position and width. The statistical uncertain- 
ties of the calibration parameters applied to MC 
are considered as a source systematic uncertainties. 
We vary each parameter by ±1 standard deviation 
to estimate the uncertainty on the final measured 
quantities. 

• MC statistics 

To determine the systematic uncertainty due to the 
limited number of MC events, we divide the MC 
samples into ten independent sub-samples and per- 
form ten pseudo-experiments. The spread of the 
unfolded ApB and measured siu^ O^g from these 
pseudo-measurements divided by vTO is assigned 
as the systematic uncertainty due to the limited 
MC statistics. 

• Electron identification efficiency 

To ensure the MC correctly models the electron 
and event selection efficiencies observed in data, 
we apply data/MC efficiency scale factors to the 
MC for forward and backward events separately. 
The bin-by-bin statistical fluctuations of these cor- 
rection factors as a function of M^g are taken into 
account and are propagated to the systematic un- 
certainties. 

• Background modeling 

To estimate the uncertainty due to the multijet 
background, we vary the reversed electron shower 
shape requirements to obtain different mass spec- 
tra of miiltijct control samples. The uncertainties 
on SM backgrounds estimated using the MC mainly 
come from the uncertainties of the energy smearing, 
data and MC efficiency scale factors, and the uncer- 
tainty of the inclusive cross section for each process. 
For the W + X inclusive background, additional 
uncertainties due to the modeling of the electron 
misidentification probability contributed by extra 



jets and the modeling of the W boson pr (obtained 
from a comparison of pythia and ALPGEN [35]) are 
also taken into account. 

• Charge misidentification 

The statistical fluctuations in the misidentification 
probability measured from data in each mass bin 
are included as a systematic uncertainty. 

The systematic uncertainties on the sin^ O^g extrac- 
tion are summarized in Table IV. The primary system- 
atic uncertainties are due to the PDFs (0.0005) and the 
EM energy calibration and resolution (0.0003). A correc- 
tion factor is introduced to account for higher order elec- 
troweak corrections which are not included in PYTHIA. It 
is determined by generating ZGRAD2 and PYTHIA samples 
and comparing the Aps distributions at the generator 
level. We find that there is a constant -1-0.0005 positive 
shift in the full sin^ O^g prediction from ZGRAD2 relative 
to the LO prediction from PYTHIA. We add this correc- 
tion factor to the extracted value of sin^ 9^,^. 

The systematic uncertainties in the unfolded ApB dis- 
tribution are listed in Table V. In addition to the com- 
mon sources listed above, uncertainties from higher order 
corrections and different SM inputs are taken into con- 
sideration. Higher order QCD, QED, and electroweak 
corrections can change the Aps predictions and thus in- 
duce additional uncertainty. We compare pythia Afb 
distributions to those of ZGRAD2 [16] with the Z/j* bo- 
son pt tuned to the RESBOS [17] prediction. RESBOS has 
the advantage of including most of the electroweak ef- 
fects with a full simulation of the non-perturbative and 
next-to-leading-logarithm (NLL) QCD ciffccts. The dif- 
ference between the two predictions is taken as a sys- 
tematic uncertainty. Different input values of sin^ in 
PYTHIA will change the kinematic and geometric accep- 
tances, and thus introduce uncertainty into the unfold- 
ing assumptions. We take sin^ — 0.232 as the default 
input value, and vary it by the measured sin^ 9^^ uncer- 
tainty (0.001). We then repeat the unfolding procedure 
and take the largest deviation from the unfolded Afb sts 
the uncertainty. 

The Qy and g\ couplings are extracted from the un- 
folded ApB distribution and thus include all of the uncer- 
tainties (statistical and systematic) that affect the ApB 
measurement. Additional uncertainties on the couplings 
from predictions with different PDF sets will be discussed 
later. 

MEASUREMENT OF sin^ Oia 

The value of sin'^ 9^^ is extracted from data by com- 
paring the background-subtracted raw Aps distribution 

with simulated Aps templates corresponding to differ- 
ent input values of sin^ 9lg. This procedure avoids the 
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TABLE IV: Uncertainties for the sin^ O^ff measurement. All uncertainties are symmetric. 



Uncertainty source 


Asin-e^ff 


Statistical 


0.00080 


Systematics 


0.00061 


PDF / Acceptance 


0.00048 


EM scale/resolution 


0.00029 


MC Statistics 


0.00020 


Electron identification 


0.00008 


Bkg. modeling 


0.00008 


Charge misidentification 


0.00004 


Higher order 


0.00008 


Total uncertainty 


0.00102 



TABLE V: Systematic uncertainties per bin for the unfolded Afb measurement. All uncertainties are symmetric. 





(GeV) EM scale/ 


Electron 


Bkg. 


MC 


PDF/ 


Charge mis- 


QCD, 


Input 


Total 






resolution 


identification 


modeling 


Statistics 


Acceptance 


identification 


QED 


sin^ eia 




50 


-60 


0.009 


0.001 


0.002 


0.011 


0.015 


0.002 


0.008 


0.002 


0.023 


60 


-70 


0.002 


0.004 


0.001 


0.014 


0.005 


0.002 


0.012 


0.001 


0.019 


70 


-75 


0.003 


0.003 


0.001 


0.015 


0.003 


0.002 


0.012 


< 0.001 


0.020 


75 


-81 


0.002 


0.001 


< 0.001 


0.008 


< 0.001 


0.001 


0.010 


0.001 


0.013 


81 


-86.5 


0.002 


< 0.001 


< 0.001 


0.005 


0.008 


< 0.001 


0.005 


0.001 


0.011 


86.5 


-89.5 


0.001 


< 0.001 


< 0.001 


0.003 


0.002 


< 0.001 


0.002 


0.001 


0.004 


89.5 


-92 


0.001 


< 0.001 


< 0.001 


0.002 


0.002 


< 0.001 


0.002 


0.001 


0.003 


92 


-97 


0.001 


< 0.001 


< 0.001 


0.002 


0.002 


< 0.001 


0.004 


0.001 


0.005 


97 


-105 


0.002 


< 0.001 


< 0.001 


0.005 


0.007 


0.001 


0.008 


0.001 


0.012 


105 


-115 


0.002 


0.001 


< 0.001 


0.008 


0.002 


0.002 


0.011 


0.001 


0.014 


115 


-130 


0.003 


0.002 


0.003 


0.010 


0.004 


0.004 


0.011 


< 0.001 


0.017 


130 


-180 


0.001 


0.005 


0.002 


0.006 


0.014 


0.011 


0.013 


0.001 


0.024 


180 


-250 


0.002 


0.002 


0.002 


0.003 


0.010 


0.011 


0.006 


< 0.001 


0.017 


250 


-500 


0.001 


0.001 


0.002 


0.008 


0.001 


0.019 


0.004 


< 0.001 


0.022 


500 


-1000 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.016 


< 0.001 


< 0.001 


0.016 



increase of the systematic uncertainty of the measure- 
ment introduced by the use of the unfolding procedure 
and maximizes the statistical significance of the final 
result. Although variations in sin^ have some ef- 
fect over the full mass range 50 < M^e < 1000 GeV, 
the central value is predominantly determined by the 
events in the Z pole region, where the statistics are high- 
est and the effects of background are smallest. Using 
events in the range 70 < Mee < 130 GeV, we measure 
sin^ = 0.2304 ± 0.0008 (stat.) ± 0.0006 (syst.) using 
PYTHIA. We then include higher order electroweak cor- 
rections using the zgrad2 program. Taking into account 
the effect of higher order corrections results in a central 
value of sin^ O'^g = 0.2309±0.0008 (stat.)±0.0006 (syst.). 
We also check the sin^ O^.^ predictions using zgrad2 and 
ZFiTTER [36] using the same input SM parameters and 
find the two results are consistent. Higher order elec- 
troweak and QCD corrections included in zfitter and 
not implemented in ZGRAD2 have a negligible impact 
on the sin^ 6'^^^ measurement. Therefore, our measured 
sin^ can be directly compared with the values mea- 
sured by the LEP and SLD Collaborations [15]. The 



comparison is shown in Fig. 6. The most precise measure- 
ments are the LEP 6-quark forward-backward asymme- 
try, A^^, the SLD left-right asymmetry, A;r(SLD), the 
LEP r-lepton polarisation measurement, Ai{Pr), and the 
SLD Icpton asymmetry, A^g. Our result is more precise 
than the LEP combined inclusive hadronic charge asym- 
metry measurement, and comparable in precision with 
the LEP c-quark forward-backward asymmetry ^4^^ . 



MEASUREMENT OF THE UNFOLDED Afb 
DISTRIBUTION 

The final unfolded Afb distribution using both CC 

and CE events is shown in Fig. 7 and Table VL together 
with PYTHIA and zgrad2 predictions. Because of the 
migration between mass bins, the correlation matrix is 
important for events near the Z pole region. The corre- 
lation coefficients are shown in Table VIL In the mass 
bins 130-180 and 250-500 GeV small deviations (< 2 
standard deviations) are observed. The x^/d.o.f between 
data and prediction is 15.3/15 for pythia, and 12.8/15 
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A,(P) 
A„ (SLD) 



I 



fb 

0, c 
fb 

had 



(D0), 5.0 fb^ 



_L 



Average 0.231 53 ± 0.0001 6 

0.23099 + 0.00053 
0.23159 + 0.00041 
0.23098 ± 0.00026 



0.23221 ± 0.00029 
0.23220 ± 0.00081 
0.2324 ±0.0012 

0.2309 + 0.0010 



_L 



_L 



_L 



_L 



J 



0.228 0.23 0.232 0.234 0.236 0.238 



FIG. 6: Comparison of measured sin^ 6^^ with results from 
other experiments. The average is a combination of A^pg, 



0,b 



aO,c 



FB-' ^'-FB^ 



and QpB measurements from 



the LEP and SLD Collaborations. 



for ZGRAD2. 



MEASUREMENT OF g'!^'-'^^ AND g^^"*^ FROM THE 
UNFOLDED DISTRIBUTION 



We extract the individual quark couplings by compar- 
ing the unfolded ApB distribution to templates generated 
with RESBOS for different values of the Z-light quark cou- 
plings. To determine g^'''' and g^'^\ the couplings of 
electrons to Z bosons are fixed to their SM values and 
sin^ is fixed to the global fit value 0.23153 [15]. A two- 
dimensional fit [38] is used to constraint the couplings, 
and a four-dimensional fit is presented as reference. The 
two-dimensional fit is performed by fixing the u quark 
[d quark) couplings to their SM values when fitting d 
quark {u quark) couplings, while the four-dimensional fit 
is performed by letting the u quark and d quark cou- 
plings vary simultaneously. The best fit values, together 
with results from other experiments, are shown in Ta- 
ble VIII. Figure 8 depicts the 68% C.L. contours of the 

fit and the contours of the theoretical uncertainty from 
the PDF uncertainties determined using the CTEQ pre- 
scription [5]. The correlation coefficients between gy, g\, 
gfr, and 5^ are shown in Table IX, without the PDF un- 
certainty included. The comparisons between different 
measurements from LEP [15], HI [37], CDF [21], and 
DO are shown in Fig. 9. Because of the high statistics 
of our data sample, and the reduced ambiguity in the 
quark content of the initial state, these are the world's 
most precise direct measurements of gy, g\, gy, and 5^ 
to date. 



0.5 







D0 5.0 fb * 1 

i 1 


ij- 


1 

- PYTHIA 




ZGRAD2 




1 Statistical uncertainty 




1 Total uncertainty 


7 , , , , 1 





50 70 100 300 500 



M^, (GeV) 

FIG. 7: Comparison between the unfolded Afb (points) and 
the PYTHIA (solid curve) and ZGRAD2 (dashed hne) predic- 
tions. The boxes and vertical lines show the statistical and 
total uncertainties, respectively. 



CONCLUSIONS 

We have measured the forward-backward charge asym- 
metry in pp — Z/7* — e+e~ events and extracted 
sin^ O^g, gy^"^^ and g^'''^ using 5.0 fb~^ of integrated lumi- 
nosity collected by the DO experiment at ~ 1.96 TeV. 
The measured forward-backward charge asymmetry in 
the range 50 < Mee < 1000 GeV agrees with the theo- 
retical predictions. The measured sin^ O'^g value can be 
directly compared with the LEP and SLD results, and the 
overall sin^ O^g uncertainty for light quarks obtained is 
smaller than the combined uncertainty in the LEP mea- 
surements of the inclusive hadronic charge asymmetry. 
We also present the most precise direct measurement to 
date oig^, g% gy^ and g%. 

Although the uncertainty of our sin^ measurement 
is still larger than that of the current world average, with 
about 10 fb^^ of integrated luminosity expected by the 
end of Tevatron Run II, a combined result of CDF and 
DO ApB measurements in both dielectron and dimuon 
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FIG. 8: [color online] The 68% C.L. contours of (a) gy and g\, and (b) gy and gj^ from a two-dimensional and a four-dimensional 
fit with statistical and systematic uncertainties. The outer regions are determined by the theoretical PDF uncertainty. The 
two-dimensional correlation contours correspond to Ax^ = 2.3 for different gA and gv parameters, as obtained from two- 
parameter (shaded regions) and four-parameter (solid and dashed curves) fits. The value 2.3 corresponds to the 68% C.L. 
region in two dimensions. In the case of four-parameter fit, the curve is a projection onto the two-dimensional plane of 
the envelope of the four-dimensional Ax^ = 4.72 surface. The crosses indicate the best two-dimensional fit values, and the 
uncertainties correspond to the one-dimensional limits. 




FIG. 9: [color online] The 68% C.L. contours of (a) gy and g^, and (b) gy and gj^ measured by DO compared with other 
experiments. The LEP and CDF Collaborations performed fits with four free parameters to determine these couplings, while 
we and the HI Collaboration performed both two and four free parameters fits. The LEP results have another solution (not 
shown) which is excluded by the HI, CDF and DO results. 
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TABLE VI: The unfolded Afb distribution compared with the theoretical predictions. The first column shows the mass ranges 
used. The second column shows the cross section weighted average of the invariant mass in each mass bin derived from pythia. 
The third and fourth columns show the Afb predictions from PYTHIA and ZGRAD2. The last column is the unfolded Afb, 
where the first uncertainty is statistical and the second is systematic. 



Mee 


(GeV) (M 


,e) (GeV) 


Predicted Afb 

PYTHIA ZGRAD2 


Unfolded Afb 


50 


-60 


54.5 


-0.293 


-0.307 


-0.305 ± 0.036 ± 0.023 


60 


-70 


64.9 


-0.426 


-0.431 


-0.431 ±0.020 ±0.019 


70 


-75 


72.6 


—0.449 


—0.452 


-0.415 ±0.015 ±0.020 


75 


-81 


78.3 


-0.354 


-0.354 


-0.343 ± 0.011 ± 0.013 


81 


-86.5 


84.4 


-0.174 


-0.166 


-0.168 ± 0.006 ± 0.011 


86.5 


-89.5 


88.4 


-0.033 


-0.031 


-0.028 ± 0.003 ± 0.004 


89.5 


-92 


90.9 


0.051 


0.052 


0.054 ± 0.003 ± 0.003 


92 


-97 


93.4 


0.127 


0.129 


0.129 ± 0.003 ±0.005 


97 


-105 


99.9 


0.289 


0.296 


0.305 ± 0.007 ± 0.012 


105 


-115 


109.1 


0.427 


0.429 


0.435 ± 0.014 ± 0.014 


115 


-130 


121.3 


0.526 


0.530 


0.576 ±0.021 ±0.017 


130 


-180 


147.9 


0.593 


0.603 


0.654 ± 0.022 ± 0.024 


180 


-250 


206.4 


0.613 


0.600 


0.623 ± 0.043 ±0.017 


250 


- 500 


310.5 


0.616 


0.615 


0.479 ± 0.068 ± 0.022 


500 


- lOOO 


."xSLl 


0.()1() 


O.Gl") 


0.812 = 0.171 = 0.016 



TABLE VII: Correlation coefficients between different Mee bins. Only half of the symmetric correlation matrix is presented. 



Mass bin 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
1") 



1 



10 11 12 13 14 15 



1.00 0.24 0.03 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
1.00 0.39 0.06 0.02 0.02 0.02 0.01 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
1.00 0.46 0.11 0.04 0.03 0.02 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
1.00 0.51 0.15 0.07 0.04 0.01 0.00 0.00 0.00 0.00 0.00 0.00 
1.00 0.72 0.32 0.11 0.01 0.00 0.00 0.00 0.00 0.00 0.00 
1.00 0.78 0.40 0.03 0.00 0.00 0.00 0.00 0.00 0.00 
1.00 0.80 0.13 0.01 0.00 0.00 0.00 0.00 0.00 
1.00 0.48 0.04 0.00 0.00 0.00 0.00 0.00 
1.00 0.38 0.03 0.00 0.00 0.00 0.00 
1.00 0.28 0.01 0.00 0.00 0.00 
1.00 0.16 0.00 0.00 0.00 
1.00 0.06 0.00 0.00 
1.00 0.06 0.00 
1.00 0.02 
1.00 



channels has the potential to substantially impact the 
world average value of sin^ 6^ff. In addition to a reduc- 
tion of the dominant statistical uncertainties, many of the 
systematic uncertainties have a strong statistical compo- 
nent or will decrease with higher statistics, for example 
the electron energy scale and rcsohition. To match the 
precision of the current world average of sin^ O^g, the the- 
oretical uncertainty due to PDFs need to be reduced in 
similar proportions as the experimental uncertainties of 
the measurement. 
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TABLE VIII: Measured Qy ' and values from different experiments compared with the SM predictions. The DO results 
are derived from best two-dimensional and four-dimensional fit, given with their total uncertainty. 





9v 






9a 


9i 


9i 



DO (2-D) 0.501 ±0.061 0.202 ± 0.025 -0.477 ±0.112 -0.377 ± 0.081 
DO (4-D) 0.501 ± 0.110 0.201 ± 0.112 -0.497 ± 0.165 -0.351 ± 0.251 
CDF[21](4-D) Q.Ultlill -0.016t°j!! -0.226lO;^« 
HI [22] (4-D) 0.56 ± 0.10 0.05 ±0.19 -0.77 ±0.37 -0.50 ± 0.37 
LEP[15](4-D) 0.471°;°! 0.241°;?? -0.52l°;°o| -0.33l°;°? 
SM [16] 0.501 0.192 -0.502 -0.347 



TABLE IX: The correlation coefficients between gy, gy and gj^. Only half of the symmetric correlation matrix is presented. 





9\ 


9"^ 


9a 


9v 


9\ 


1.000 








9v 


0.470 


1.000 






9i 


0.201 


-0.606 


1.000 




9'v 


0.217 


0.925 


-0.813 


1.000 
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